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Stochastic perturbation of two-level atoms strongly driven by a coherent light field is analyzed 
by the quantum trajectory method. A new method is developed for calculating the resonance 
fiuorescence spectra from numerical simulations. It is shown that in the case of dominant incoherent 
perturbation, the stochastic noise can unexpectedly create phase correlation between the neighboring 
atomic dressed states. This phase correlation is responsible for quantum interference between the 
related transitions resulting in anomalous modifications of the resonance fiuorescence spectra. 
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I. INTRODUCTION 

Quantum interference is one of the most intriguing phe- 
nomena of quantum mechanics. Over the past decade 
several effects in atom-hght interaction which have their 
origin in quantum interference have been predicted and 
demonstrated experimentally ||l||. Some characteristic 
examples are reduction and cancellation of absorption 
lH, ^, ^, ^ and spontaneous emission 1^, ^, ^, |l^ , and 
narrow resonances in fluorescence p^f T^ A prerequisite 
of quantum interference between the transition channels 
is the existence of some stable time correlation of the 
atomic system under consideration. A possible way of 
achieving such correlation is the application of coherent 
coupHng in a multi-level atomic system. Although some 
interference effects have also been found in two-level sys- 
tems interacting with two light beams ||l^, so far quan- 
tum interference has been observed exclusively in at least 
three- level systems. 

Generally, various incoherent perturbations destroy 
the phase correlation between the states involved in the 
interfering transition pathways and the coherently in- 
duced quantum interference disappears. However, under 
special circumstances, even incoherent perturbation can 
be responsible for quantum interference. For example, 
quantum interference can appear in three-level systems 
due to collisions. Such effects are known as pressure- 
induced extra resonances and have been studied in four- 
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wave mixing ||l 4 |l5| . 

Recently, in an experiment with coherently driven two- 
level atoms, anomalous resonance fluorescence spectra 
were found when the collisional relaxation rate exceeded 
the Rabi frequency |[l6). The spectra were of the form 
of a pressure-broadened line with a narrow, not colli- 
sionally broadened dip. These results, unexpected in a 
collisionally perturbed two-level system, were interpreted 
as a consequence of quantum interference between differ- 
ent dressed-state transition channels. In Ref . ||l6|| , it was 
pointed out that these effects can also occur in the case 
of a non-monochromatic, e.g. phase-diffusing, laser fleld. 
Indeed, resonance fluorescence spectra with the phase- 
diffusing laser fleld have been calculated by Peng Zhou 
et al. [13 who obtained the same results as GawHk et 
al. What in both cases appears essential for obser- 

vation of quantum interference and anomalous spectra 
is that the incoherent perturbation (collisions or phase 
diffusion of the light fleld) dominates over the Rabi os- 
cillations. 

These two examples raise an intriguing question how 
can a stochastic noise lead to stable time correlation re- 
sulting in quantum interference in two-level systems. 

Quantum trajectory methods are widely used powerful 
tools for treating the stochastic evolution of open quan- 
tum systems |jl§ |l9[ |o[ |2^, |2|, H]. They can provide 
the solution of any master equation that is of Lindblad 
form p^ , |2^ , p^ , p7| |. Moreover, individual quantum tra- 
jectories, as state evolutions conditioned on particular 
sequence of observed events, make it possible to reveal 
phase correlations in the given system. 

While resonance fluorescence spectra can be ade- 
quately modeled by the master equation, we flnd the use 
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of quantum trajectories provides more physical insight. 
Besides, this method allows to study the stochastic evo- 
lution of the atomic system and, eventually, reveal the 
phase correlation of its dressed states. For that reason, 
in this paper we analyze in detail the system of stochasti- 
cally perturbed two-level atoms applying both the master 
equation and the quantum trajectory method to explain 
quantum interference effects and the underlying physical 
processes. 

The paper is organized as follows: In Section II we in- 
troduce our model for the system of a coherently driven 
and incoherently perturbed two-level atom. In Section 
III the method of quantum trajectories is applied to the 
system and a new method is developed for calculating 
the spectrum from the simulation results. Section IV 
presents the numerical results for the spectra and com- 
pares them with the analytical solution of the master 
equation. In Section the phase difference between the 
dressed states of an atom is analyzed and the phase corre- 
lation is revealed. It is shown that the phase correlation 
emerges as the noise magnitude increases and the related 
quantum interference effect is discussed. 



II. THE MODEL 

The system of incoherently perturbed and coherently 
driven two-level atoms can be modeled in several ways. 
Here we make a rather general assumption that the 
stochastic perturbation is responsible for fluctuations of 
the atomic resonance frequency which obeys the Gaus- 
sian statistics. In particular, such fluctuations may result 
from e.g. elastic, dephasing collisions. 

In our model the Hamiltonian of the strongly, coher- 
ently driven atom subjected to stochastic perturbation 
has the form 



H 
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(1) 



in the interaction picture, where uja{t) = cJq + 5uja{t) 
is the fluctuating atomic transition frequency, t^L is the 
frequency of the laser, Q. the Rabi frequency, and , , 
S~ are the atomic operators deflned in the excited state 
(|e)) - ground state (l^)) basis: 



qz_^ { 1 0\ + /01\ /OO 

"2I -ij'-^ ~Uoj'^ "UO 
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We assume that the noise in the transition frequency sat- 
isfles 



(3) 



with fluctuating phase, if the phase drift is neglected |^ . 
In such case, F represents the laser linewidth. 

The time-evolution of the system deflned by the Hamil- 
tonian in Eq. (|l]) is described by the following master 
equation, taking into account also the spontaneous emis- 
sion processes: 



where 
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and (Hal) is the mean atomic Hamiltonian obtained by 
averaging over the stochastic noise of Eq. (H), and 7 is 
the natural linewidth of the atom. 

It is natural to introduce the dressed-state basis in 
which the atomic Hamiltonian Hal is diagonal: 



|1) 
|2) 



cos0|g) + sin9|e) 
sin6|g) + cos6|e) 



(6a) 
(6b) 



where 
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and A — uja — ujl is the laser detuning. 

The mean atomic Hamiltonian {Hal) in the dressed- 
state basis |1), 1 2) can be written as 



where 



{Hal)=Ei\1){1\+E2\2){2\, 



(7) 



The effect of incoherent perturbation on the pur e states 
can be determined from the Lindblad form ( |5c|) of the 
master equation in Eq. (Q). The action of the oper- 
ator 2\/TS^ corresponds to an event generated by the 
stochastic noise. Without detuning, this operator gener- 
ates transitions between the dressed states |1) and |2): 



2VTS'\1) = \/f|2), 

2Vrs'\2) = Vt\i). 



(8a) 
(8b) 



Another way of observing the effect of the stochastic 
noise is to transform the time-dependent Hamiltonian in 
the Langevin equation into the dressed-state basis. This 
must be done carefully since parameter Q in the deflni- 
tion of the dressed states becomes time-dependent in this 



where F stands for the magnitude of the stochastic noise. 
If the noise is due to collisions, this quantity is the col- 
lision rate between the atoms. This model can also de- 
scribe the system of two level atoms driven by a laser fleld 
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to the first order in Suja{t). Hamiltonian Hal is diagonal 
in the time-dependent dressed-state basis, thus 



Hal = Ei\l,t){l,t\+E2\2,t){2,t\ « 
^E,\l){l\+E2\2){2\ + 

+S^ait)l^f^i-\1){2\-\2){1\) 
'Su;ait)l^r^^i\l){2\ + \2)(1\) ^ 



(10) 



2 A2 + f}2 

= {Hal) -h5uja{t) 



:(|1>(2| + |2)(1|). 



This also shows that the stochastic noise generates tran- 
sitions between the dressed states |1) and |2). 

The master equation has the following form in the 
dressed-state basis: 



_d 
dt 
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where = pn - P22, T' = T - 7/4, and pn, pi2, P21, 
P22 are the matrix elements of the density operator in 
the dressed-state basis. The matrix element p2i is the 
complex conjugate of pi2, as p is Hermitian. For reso- 
nant excitation (A = 0) the stochastic noise couples the 
|1) and 1 2) states and increases the relaxation rate of the 
spin z component. In this case the dressed states become 
independent of the Rab i frequency, and Eq. ( |lla| ) is un- 
coupled from Eq. ( jllbl ). This is however not the case in 
the general, non-resonant case. 



III. NUMERICAL SIMULATION 

In the system of coherently driven stochastically per- 
turbed two- level atoms, quantum interference eflects can 
be seen in the resonance fluorescence spectra ||l^, 0. 
The resonance fluorescence can be described by transi- 
tions between appropriate dressed states of the atom. If 
spectral modifications are due to quantum interference, 
some time correlation should exist between the dressed 
states of the atom involved in the interfering transition 
channels. 

For analyzing time correlations in a quantum system, 
quantum trajectory methods are particularly appropri- 
ate. These methods are based on the simulation of quan- 
tum trajectories, that are individual realizations of the 



evolution of the system conditioned on particular se- 
quence of observed events. By tracking the time evolu- 
tion of a single quantum trajectory, the time correlations 
can be revealed. 

We apply the quantum trajectory method of Ref. ||l^ 
for simulating the time evolution of the coherently driven, 
stochastically perturbed two-level atom. In this system, 
a single quantum trajectory evolves coherently according 
to the Hamiltonian of Eq. (^ , interrupted by incoherent 
gedanken measurements due to noise events and sponta- 
neous emission. The evolution of the density operator of 
the system is obtained by averaging the density opera- 
tors of the individual quantum trajectories. The resulting 
density operator is the solution of the master equation of 
Eq. (|). 

The accuracy of the simulation is limited by two fac- 
tors: the length At of the time step and the number N of 
the simulated quantum trajectories. At should be much 
less than the characteristic time of any process in the 
system. N should be large enough to obtain the right 
ensemble averages for the density operator at the given 
stochastic noise magnitude. In our simulations N was 
approximately 5 • 10^. 

Within dipole approximation the resonance fiuores- 
cence spectrum S{uj) can be calculated as the real part 
of the two-time correlation function 



POO 

rf(w)=lim / exp{^iujT){S+it + T)S-{t))dT, 



for an arbitrary initial condition: 

S'(w) =Rerf (w), 



(12) 



(13) 



where uj is the detuning of the emitted light from ujl- 
There are different methods in the literature for obtaining 
the spectrum using a numerical simulation | p^ , p5| , p6[ . 
One kind of them simulates not only the atom but also 
the quantized electromagnetic field ||2^. Such methods 
seem to be excessive when the field can be treated clas- 
sically. 

The method presented by Dalibard et al. |Q simulates 
only the atom, and obtains the spectrum by calculat- 
ing two-time averages and taking their Fourier-transform. 
The computation time of this method increases as 1/At'^, 
where At is the time step of the simulation, because for 
each time step an additional simulation is started to cal- 
culate the two-time averages. That can be time consum- 
ing in the case when large number of quantum trajecto- 
ries are simulated and small time steps are used. This is 
the situation in our problem when we simulate the system 
in the high noise magnitude regime. 

The question arises whether it is possible to develop a 
method which simulates only the atom, without the need 
of starting extra simulations for calculating the two-time 
averages. Below, we briefiy outline the essentials of our 
novel method for spectrum calculation. More details will 
be published elsewhere. 
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Let us consider a general two level atom-field system. 
Let p{t) be the density operator of the whole system, A 
an operator in the Schrodinger picture acting only on 
the atom, and U{t) the unitary time-evolution operator. 
Then 

Tr{Ap{t)) = Tv[AU[t)p{m\t)) = 

= Tv{U\t)AU{t)pm^Tv{A{t)pm, 

where A{t) is the operator A in the Heisenberg picture. 
One can define a time-dependent A'{t) operator for which 

Tr{Ap{t)) - TYA{ApA{t)) = TrAiA'it)pAm, (14) 

where pA{t) = TrLp(i) is the reduced density operator 
of the atom. A'{t) depends also on the laser field. Let 
i?i(0) be a set of reduced density operators of the atom 
that form a C-linear basis in the set of the operators 
acting on the atom. In the case of a two-level atom the 
basis consists of four elements. These basis elements may 
evolve also, their value at time t is denoted by Ri{t). 
Using this basis, any operator X which acts on the atom 
can be written in the form 



i 

The coefficients can be expressed as 



(15) 



(16) 



where T^j = Tr(i?i(0)i?j (0)). Matrix T is invertable since 
the operators Ri (0) form a basis and they are linearly in- 
dependent. Using Eqs. (p^, (|l^ ) and ( |l6|) the following 
form can be derived for the A'{t) operator: 

A'{t) = ^(r-i),feTr(i?,(i)i?.(0))A.(0)i?,(0) (17) 



where (0) is expressed by the atomic operator A using 
Eq. @ as 

A.(0) = (T-i).,Tr(Ai?,(0)). 

Equation ( |l4|) h olds for all density operators p(0) with 
A'{i) ofEq^): 

Tv{A{t)pm=TrA{A'{t)pAm. (18) 

Having an operator B{t) in the Heisenberg picture such 
that B{Q) acts only on the atom, B{Q)p{Q) can be ex- 
pressed as a C-linear combination of density operators. 
The linearity of the trace in Eq. (tt8) yields 



Tr (A(i)S(0)p(0)) = Tr {A' {t)B{Q)pAm • (19) 

For calculating two-time correlation functions of the form 
Tt A{t)B{t')p{Q), the above equation can be modified by 
using the cyclic property of the trace: 

TTA{t)B{t')p{{)) = i:rA{t)U{-t')B{{))U{t')U{^t') ■ 

■p{t')U{t') = TYU{t')A{t)U{-t')B{Q)p{t') = 
= TrA(t-t')S(0)p(t') = 

= i:TAA'{t-t')B{0)pA{t'), (20) 



where pAit') is the reduced density operator in the 
Schrodinger picture at time t' . 

Let us apply the general expressions presented above 
to the atomic operators S'^ and S~ . Calculating the 
correlation function, the quantity Tr S^{t -I- r)S'~(i)p(0) 
should be determined from the simulation. Using 
Eq. @, 

lim TYS+{t + T)S-{t)p(Q) = lim Tr5+(T) • 

t—^oo t—*oo 

■S-{Q)p{t) = Tr^ S+\t)S-{Q)pa{^). (21) 

In order to obtain operator 5'+'(r) from the simulation 
one needs to choose a basis consisting of density opera- 
tors, according to Eq. ( |l5|) and start independent simu- 
lations using the elements of this basis as initial states. 
In our simulation we choose the density operators that 
in the dressed-state basis |2) are: 



i?i(0) 



1 




^2(0) = - 



1 1 
1 1 



^3(0) = - 



1 

1 -1 
-1 1 



(22) 



For all time steps of the simulation we calculate the 
S'"*"'(t) operator using Eq. ( [l7|) and record it for later 
use. After the simulation has been completed, i.e. the 
time t has reached its final value T, the correlation func- 
tion defined by Eq. (^2|) is calculated numerically by 
evaluating the expression 

T 

rf (w) = ^exp(-it^T)Tr5+'(T)5-(0)^^At, (23) 

where the su mma tion is done over all time steps between 
and T and p[T) is the average of all four Ri{T) density 
operators. 

The advantage of this method is that there is no need 
to start a new simulation in each time step, and it is suf- 
ficient to simulate only the atomic system for obtaining 
the spectra. 



IV. THE FLUORESCENCE SPECTRUM 

In order to check our numerical results, we compare 
them with the spectra calculated analytically. After de- 
termining the time evolution of the averages of the Block 
vector components {S^{t)), (5+(t)), (S^{t)), the quan- 
tum regression theorem is used for expressing the two- 
time average {S+{t + r]S'"(t)) in Eq. (|l|) as a function 



of one-time averages 
following: 



28|. The Bloch equations are the 



{S+it)) = 



-inns- it)) 



1 



After some calculation we obtain 
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injiu - iA + 2r + ^7) {-^K3 - K1K3) 

{iu + 7) {{iuj + 2T+ ^7)2 + A2) + n-^{iuj + 2r + ^7) ^ 

(^^2 ^ (^^ ^ ^)(^^ - iA + 2r + ^7)) + j^l - j^2J^3 

{iuj + 7) ((iw + 2r + i7)2 + A2) + n^{ioj + 2r + 17) 



ir!2 (-i^gi^s) 



where 



K2 
K3 



i7((2r+i7)' + A2 



7 ((2r + ^7)2 + A2) + r!2(2r + 17) ' 



^jin (iA - 2T 



7 



7 ((2r + ^7)2 + A2) + r!2(2r + ^7) ' 

^-fifi {iA + 2r + ^7) 
7 ((2r + i7)2 + A2) + r!2(2r + 17) • 



For the special case of no detuning (A = 0) the correlation function has the form 

1 



(4r + 7) 72)172 



- 2i{T + 17)) (iia2(2w2 „ 4j^^r" - a2)) 
+uj^{2a^ - 2^22(^2 + 72)) + icj ((37 + 4r")r2 



217^ 



4r"a'* - 27(2172 - 72)172 



217'' 



(24) 



(25a) 
(25b) 
(25c) 

(26) 



where 



a2 ^72 + 4r7 + i72, r" = r' + 7, r' = r--7. 



(27) 



The correlation function of Eq. ( ^61 ) can be split into 
the sum of three functions: 



rf(c.) 



A, 



A- , Ao 



U! — OJ — S- OJ — Sq 



where 



s± = ir"±iVr'2-i72, 

So = 2ir' + ij, 
A, ..172 + 



(28) 



(29a) 
(29b) 



2ia2(s_|_ — s_)(zs_ + 7) 
.o2 -47'r"(r+ I7) ±Z7a2(s+ - S-) 



2i(s+ — s_)a*(is_ + 7) 



Ao = -i 



172 + (4r + 7)7 

2a2 



If r' < 17, the spectrum has the form 

Aoso RcA+iv - Rc{A+s*^) 



'^ + ko|2 ' (cj + v/172 - r'2) + r" 

ReA^uj - Re{A^s*_) 

{lu - Vf^2 _r'2) + r"2 ' 



(30) 



showing that the centers of the two latter Lorentzians are 
displaced by zt\/ 172 _ j-gjg^^jyg ^jjg j^ser frequency. 
Together with the first Lorentzian at the laser frequency 
they form the Mollow triplet 

In the other case when F' > 17, all the Lorentzians 
are centered at zero frequency, corresponding to in 
the Schrodinger picture, but one of them has a negative 
coefficient, resulting in a dip in the spectrum. 

For the F' > 17 case, one obtains the following expres- 
sion for the spectrum: 



a;2 + |s+|2 



AqSq 



\so\ 



(31) 



IS 



where A+s+ and AoSo are always positive and A^s 
negative. 

In the following we show our numerical results together 
with the analytical spectra. Fig. ^ presents the reso- 
nance fiuorescence spectrum of the atom irradiated by 
a resonant (A = 0), strong (17 ^ 7) laser field with low 
noise (17 > F). The spectrum exhibits a three-peak struc- 
ture, but with a suppressed and broadened central peak 
compared to the standard Mollow triplet. As the noise 
increases, the central peak disappears and for F nearly 
equal to 17 we get a two-peak structure with a relatively 
broad dip, as depicted in Fig. ||. When the noise mag- 
nitude is much larger than the Rabi frequency, the dip 
becomes very narrow, as shown in Fig. g[ The width of 
the dip is pr oportional to the value of the parameter |s_| 
of Eq. (29a), which approaches the natural Hnewidth 7 
when F' ^ 17, as shown in Fig. |4[ For large detuning 
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(A ^ ft) and low noise (F ^ f2), a two-peak spectrum is 
obtained with an asymmetric Fano-like structure at the 
center, as depicted in Fig. ||. Increasing the noise mag- 
nitude the Fano-hke peak transforms to an asymmetric 
Fano profile, a narrow dip on the side of the broadened 
part of the spectrum due to stochastic noise, and a nar- 
row peak on the other side next to the dip (Fig. ^). 

The analytical and numerical results are in very good 
agreement as one can see in the presented figures which 
positively verifies correctness of our numerical simula- 
tions. 



V. PHASE CORRELATION AND QUANTUM 
INTERFERENCE 

The narrow dip in the spectrum (in the case of high 
noise magnitude, resonant excitation) and the asymmet- 
ric Fano profile (in the case of large detuning and high 
noise magnitude) are signatures of quantum interference 
in the stochastically perturbed system. The quantum 
interference emerges if long-time phase correlation exists 
between quantum states connected by different transition 
channels. As seen in section II, the stochastic noise gen- 
erates transitions between the dressed states. It would be 
interesting to check whether this couphng has any effect 
on phase correlation between them. 

The phase difference A(j> between the dressed states is 
defined as follows: 



1$) =aie^*Hl)+a2e"^l2), : 



(32) 



where |<i>) is a pure state of the atom, while |1) and 
1 2) are the dressed states defined in Eqs. ( |6al ) and (|6b|). 
The phase difference can be calculated straightforwardly 
from a single quantum trajectory. It is found that the 
phase difference behaves differently in the low and high 
noise magnitude regimes. In Fig. 0a the noise magni- 
tude is much less than the Rabi frequency, Rabi oscil- 
lations are rarely disrupted by noise events, hence the 
phase difference is essentially hnearly dependent on time: 
A(j){t) = 2flt. Consequently, the shape of the phase dif- 
ference as the function of time shows no structure. When 
the noise magnitude increases, as depicted in Fig. 0b, the 
uniform shape changes to a picture showing some struc- 
ture of gaps appearing from time to time between and tt 
values of the phase difference. For high noise magnitude 
(Fig. 0c), the phase difference tends to stabihze around 
values and tt for some time intervals. 

In order to characterize the observed phenomena quan- 
titatively, we introduce the correlation function of cos Acj) 
by the definition 



Ccos(t) 



c / (cos A0(i + r) — cos A0) 
/t=o 



X (cos A(j){t) — cos Acj)) dt, 



(33) 



where cos Acj) is the mean value of the cosine of the phase 
difference for the simulated time interval and c is a nor- 



malization constant fixed by the condition Ccos(O) = 1. 
The correlation function of sin Acf>{t) is defined similarly. 

The Ccos(''') function is shown in Fig. Js] for the same 
parameter values as those used in Figs. ma^c. The quali- 
tative picture of emerging correlations as the noise mag- 
nitude increases is now backed up by the widening of the 
correlation functions. On the other hand, the correla- 
tion of the sine of the phase difference, Csin(T) (defined 
similarly as Ccos(''')) tends towards a S-\ike shape when 
the noise increases (Fig. ||), so when T strongly exceeds 
fl, smAcj>{t) remains uncorrelated. This means that the 
phase difference is locked to values and tt for some time 
intervals, though it spans a phase interval no less than j 
around these phase values. 

It is interesting that the widths of the correlation func- 
tions are related with the observed narrow dips in the 
spectra. The full width at half of the maximum value 
(FWHM) is a good measu re o f the widths of Ccos('''), 
and \s-\ as defined in Eq. (29a) describes well the spec- 
tral dip width. These two quantities are presented in Fig. 

The FWHM of Ccos(r) is roughly proportional to the 
reciprocal of the width of the dip in the spectrum across a 
wide range of parameter sets, so the observed phase cor- 
relation is indeed responsible for the narrowness of the 
dip in the spectrum. 

The stabihzation of the phase difference between the 
dressed states of the stochastically perturbed and co- 
herently driven two-level atom is the underlying physi- 
cal process which makes the quantum interference pos- 
sible. This stabilization supports the following inter- 
pretation first suggested for collisional and phase noise- 
induced quantum interference effects in resonance fiuo- 
rescence spectrum in Ref. jl^. Resonance fiuorescence 
of a strongly-driven two-level atom is emitted in cas- 
cade transitions downward the ladder of the dressed-state 
doublets. Fig. |ll| shows two adjacent doublets and all 
possible spontaneous and noise-induced transitions be- 
tween the dressed-atom states. According to Eqs. (||) 
and (|l0|), noise events generate transitions between the 
dressed states |1) and 12) and couple them as indicated by 
double arrows in Fig. As we have seen in the previous 
section, in the noise-dominated regime, i.e. when F > 17, 
the phase difference between dressed doublets tends to 
stabilize for some time intervals due to frequent noise 
events. Moreover, according to Eq. (|3T|), the resonance 
frequencies of all fiuorescence contributions are the same 
in this regime. Among several possible emission channels 
there are two pairs: \l,n) —f |2, n — 1) 1) and 

\2,n) ^ (or \2,n) ^ \l,n) |2,n-l) 

and |2,n) 1) — > |2,n— 1)) that differ exclusively 

by time ordering between collisional mixing and photon 
emissions. Photons emitted along these channels are in- 
distinguishable, so their interference is possible. Due to 
opposite signs of the relevant matrix elements this inter- 
ference is destructive and creates a dip in the line center. 
On the other hand, other emission channels are not that 
equivalent, hence the corresponding photons cannot in- 
terfere and contribute to non-zero intensity at uj = 0. 
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This interference is similar to that seen by Schrama et 
al. in photon correlations of the well resolved Mol- 
low triplet components in the opposite limit when F < fi. 

VI. CONCLUSION 

We have applied the quantum trajectory method to 
the system of two-level atoms strongly driven by a co- 
herent light field and perturbed by stochastic noise. We 
have developed a new method for obtaining the reso- 
nance fluorescence spectra from numerical simulations. 
This method is especially advantageous for physical sys- 
tems where the noise dominates the dynamics, and one 
needs to simulate many quantum trajectories using small 
time steps. The simulation of a single quantum trajec- 
tory revealed that for high noise magnitude the phase 
difference between the dressed states tends to stabiHze 



around fixed values. When calculating the resonance fiu- 
orescence spectra, narrow resonances as central dip and 
dispersive Fano-like profile occurred in the regime where 
the noise dominated the Rabi oscillations. These modi- 
fications of the resonance fiuorescence spectra are asso- 
ciated with the stabilization of the dressed-state phases 
and stochastically-induced quantum interference between 
various emission channels. 
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FIG. 1: The resonance fluorescence spectrum for low noise 
magnitude and strong laser fleld (F/fl = 0.2, 7/f2 = 0.05) in 
the case of no detuning (A = 0). 
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FIG. 2: The resonance fluorescence spectrum for noise mag- 
nitude comparable to Rabi frequency and strong laser fleld 
(F/n = 1.1, 7/a = 0.05) in the case of no detuning (A = 0). 
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FIG. 3: The resonance fluorescence spectrum for high noise 
magnitude and strong laser fleld (F/f2 = 6, = 0.05) in 
the case of no detuning (A = 0), showing a narrow dip at the 
center of the spectrum. 
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FIG. 4: The width of the dip in the case of no detuning, i.e. 
the value of shown for j/fl = 0.05. When V increases, 
the width of the dip approaches the natural linewidth 7. 




-15 -10 -5 5 10 

co/a 



FIG. 6: The resonance fluorescence spectrum for high noise 
magnitude, strong laser field and large detuning (F/fi = 3, 
7/n = 0.05, A/f2 = 3), showing an asymmetric Fano profile 
at the driving frequency. 
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FIG. 5: The resonance fluorescence spectrum for low noise 
magnitude, strong laser field and large detuning (F/fi = 0.2, 
7/f2 = 0.05, A/f2 = 3), showing a Fano-like structure at the 
driving frequency. 




FIG. 7: The phase difference between the dressed states 
in the case of resonant excitation (A = 0), a) for a low noise 

magnitude (F/fl = 0.2, 7/n = 0.05); b) for a noise magnitude 
comparable to the Rabi frequency, (T/il = 1.1, 7/51 = 0.05); 
c) for a high noise magnitude, {T/Q = 5, 7/O = 0.05). The 
initial state was the excited state |e) in these simulations. 
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FIG. 8: Correlation function of cosA(j){t) for low, medium 
and high noise magnitude, in the case of no detuning and 
strong laser field [T/Q 6 {0.2, 1.1, 5}, A = 0, -y/Q = 0.05). 




FIG. 9: Correlation function of sin A0(f) for low, medium and 
high noise magnitude, in the case of no detuning and strong 
laser field [T/Q, € {0.2, 1.1, 5}, A = 0, 7/n = 0.05). The 
figure shows that sinA0(t) becomes uncorrelated for higher 
noise magnitude. 
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FIG. 10: The FWHM of the correlation function Ccos(t) is 
plotted against the reciprocal of the analytically calculated 
dip width |s- 1. 
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FIG. 11: Different dressed-state transition channels for a co- 
herently driven and stochastically perturbed two-level atom. 



